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in : ABSTRACT 

p/] ' Core-accretion planet formation begins in protoplanetary disks with the 

^ ■ growth of small, ISM dust grains into larger particles. The progress of grain 

'^ ■ growth, which can be quantified using lO/xm silicate spectroscopy, has broad 

I ' implications for the final products of planet formation. Previous studies have 

H ' attempted to correlate stellar and disk properties with the lO/xm silicate feature 

52 ■ in an effort to determine which stars are efficient at grain growth. Thus far there 

does not appear to be a dominant correlated parameter. In this paper, we use 
spatially resolved adaptive optics spectroscopy of 9 T Tauri binaries as tight as 
0.25" to determine if basic properties shared between binary stars, such as age, 
^ ' composition, and formation history, have an effect on dust grain evolution. We 

fr^ ■ find with 90-95% confidence that the silicate feature equivalent widths of bina- 

[■^ ■ ries are more similar than those of randomly paired single stars, implying that 

shared properties do play an important role in dust grain evolution. At lower sta- 
tistical significance, we find with 82% confidence that the secondary has a more 
prominent silicate emission feature (i.e., smaller grains) than the primary. If 
^ ■ confirmed by larger surveys, this would imply that spectral type and/or binarity 

j^ ■ are important factors in dust grain evolution. 



Introduction 



The core-accretion model for planet formation (e.g., iLissauer &: Stevenson! 120071 ) is a 



multi-step process beginning with the agglomeration of small interstellar medium (ISM) 



-'^The observations reported here were obtained at the MMT Observatory, a facihty operated jointly by 
the Smithsonian Institution and the University of Arizona. 
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dust grains into larger particles. Eventually, these particles reach the size of planetesimals 
and gravitationally attract each other and their surrounding gas. Planet formation must 
occur quickly, as several physical processes are able to disperse the gas and dust of the 
protoplanetary disk. Consequently, the timescale over which each aspect of core-accretion 
occurs is critical to the final outcome of the planet formation process. The first step of core- 
accretion, where small ISM dust grains agglomerate into larger particles, is of particular 
importance as it sets the conditions for all of the other core-accretion steps. 

Dust grain pro perties can be studied at a variety of wavelengths and locations in the disk 
( iNatta et al.l 120071 ). Visible and near-infrared scattered light imaging can be used to study 
the grain sizes of submicron and micron surface dust at large radial distances. Mid-infrared 
spectroscopy of the lO/xm and 20/im silicate emission features can be used to determine the 
size and composition of submicron and micron sized surface dust at small radial distances. 
Far-infrared and millimeter continuum observations probe deep into the disk mid-plane over 
the full radius of the disk and can be used to determine the size distribution and total mass 
of mm-cm dust grains. Although these observations probe different grain sizes at different 
locations in the disk, it appears that they are somewhat correlated. In particular, the equiv- 
alent width of the 10/xm silicate feature is correlated with the slope of the submillimeter 
spectral energy distribution (SEP), implying tha t dust grain growth proceeds concurrently 
through its various phases (JLommen et al.ll2010l ). Additionally, the shape of the 10/im sil- 
icate feature is correlated with its amplitude, like ly due to the sin i ultan e ous development 
of large amorphous gr ains and crystalline grains ( 
2OO3I : lApai et allboosl ). 



van Boekel et al.l l2003l : iPrzygodda et al. 



Our current understanding of dust grain growth is that it begins during t he Class II 
phase of star /planet formation, once the dust is collected into a cir cumstellar disk (JBeckwith fc Sargent 
I991I : lOssenkopf fc Henning]|l994l : lKruegel fc Siebenmorgerull994l ). At this point, dust grains 
are expect ed to coagulate and settle rapidly (~ 10^ ~ 10^ years depending on model as 
sumptions; 



DuUemond &: Dominik 2005 



Brauer et al.ll2008l ). However, Class II objects over 
0.5-10 Gyr) have a diverse se t of dust grain properties, implying 



a wide span of ages ^ „ ^ 

that dust evolution is not simply a function of age (IMeeus et a. 



2003 



Oliveira et al 



van Boekel et al.ll2005l : lApai et al.ll2005l : iKessler-Silacci et al. 



2001 



2006 : 



op 



Przvgodda et al. 



Furlan et al. 



2009 



^^^^^M ■ ± ■ ■ ■; ■ III; ■ III; 

201 ll ). Even within individual clusters, which are approximately coeval, there 



is no apparent correlation betw 'een mid-infrared dust properties and stellar ma ss, luminosity, 
stellar accretion or disk mass (jSicilia-Aguilar et al.l 120071 : IWatson et al.l |2009| ) . Over a large 
range of masses, taken between stars in different clusters with differen t ages, there appears 



to be a weak correlation b etween mass and silicate feature strength (JKessler-Silacci et al. 



20061 : iPascucci et al.ll2009l ). which might be the result of another spe ctral-type dependent 



factor, such as X-rays or the location of the silicate emission zone (JGlauser et al.l 12009 
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Kessler-Silacci et al.ll2007| ). However, the general diversity of dust properties between simi- 
lar stars of similar ages suggests that additional properties are needed to explain the evolution 
of the dust. 

One way to isolate which stellar properties are important for dust grain evolution is to 
observe the dust grain properties of coeval young binaries. Binaries share certain properties, 
such as age, composition and formation history, which are difficult to ascertain in individual 
stars. If these properties are important for dust grain evolution, then we expect the dust 
to be similar in binary pairs. However, if these properties are unimportant, then the binary 
stars will have a random distribution of dust characteristics, similar to those observed in 
single stars. Additionally, if the shared properties are important, they might be hiding 
correlations between other stellar parameters and dust grain characteristics, which can be 
uncovered by taking advantage of the coevality of binaries. 



In this paper (and including the results of ISkemer et al.ll2010l ). we present spatially re- 
solved 10/im silicate spectroscopy of 8 Taurus- AurigacI binaries, which triples the sample size 
of spatially /spectrally resolved binaries in a single cluster at mid-infrared wavelengths. Small 
diameter space-based telescopes are unable to perform high spatial resolution (< 2 — 3") ob- 
servations in the mid-infrared (>8/im) due to their large diffraction limits (additionally, the 
Infrared Spectrograph, IRS, on Spitzer is no longer operational, so mid- infrared spectroscopy 
is not currently accessible from space). Ground-based observatories have much larger diam- 
eters but suffer from high sky background, variable transmission and seeing, which explains 
why constructing even a small sample has been difficult. We describe our observations and 
reductions in Section 2 and our corrections for extinction and equivalent width measure- 
ments in Section 3. In Section 4, we perform statistical tests in an attempt to determine 
which properties affect dust grain-growth. In particular, we ask the questions (1) Is grain 
growth correlated between binary pairs? and (2) Does removing this effect reveal correlations 
between grain growth and other properties? We discuss the implications of our results and 
our conclusions in Section 5. 



2. Observations and Reductions 

We observed 7 T Tauri binaries with MMTA0/BLINC-MIRAC4 over several nights 
in 2009 and 2010 (the object names, dates, and exposure times can be found in Table 
1). All of the objects were observed with the 6.5 meter MMT Telescope and its unique 



^Taurus- Auriga is a low-mass star forming; r egion with an age of '^1-2 Myr and a distance of '^140 pc 



(JKenvon fc Hartmannlll995l : lKenvon et aLlll994) . 
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deformable-secon dary adaptive optics system (MMTAO: lLloyd-Hartll200Cll : IWildi et al.ll2003 
Brusa et al.ll2004l ). which pr oduces stable, high -Strehl (up to ~98%) point-spread-functions 



(PSFs) in the mid-infrared ( jClose et al.ll2003l ). We used the combination mid-IR nulhng- 
interferometer /imager/spectrograph, BLINC-MIRAC4 (Bracewell Infrared Nulhng Cryostat; 
Mid- I R Array Camera , Gen. 4), in its irn ager /spectrograph configuration ( iHoffmann et al. 
19981 : iHinz et al.ll200ol : Iskemer et al1l2009h . 



For each binary, we observed a bright spectroscopic standard, that doubled as a PSF 
star, before and after the binary observation. We took acquisition images at 8.7/im to put 
the standard or binary in MIRAC4's 1" xl5" slit. We then inserted the KRS-5 grism, which 
is in a filter wheel, near a pupil plane. For the binaries, we azimuthally aligned both stars 
with the slit, either based on our 8.7/im images or on optical acquisition images for the 
AO-system's wavefront sensor. The telescope's Cassegrain de-rotator was turned on to keep 
both stars in the slit. We chopped perpendicularly to the slit with an 8" throw, and nodded 
along the slit (so that the binary was visible in both nods) with ~3-5" nods, depending on 
the binary separation. 

The adaptive optics (AO) system is able to run at full-speed (550 Hz loop-speed) for 
stars brighter than V~12 mag, and sometimes fainter in good conditions. Since some of our 
targets were fainter than this cutoff, we ran the AO system slower (100-150 Hz loop-speed) on 
these binaries, and their PSF star/spectroscopic standards. Our PSF/spectroscopic standard 
stars were brighter at visible wavelengths than our target binaries, so we used neutral density 
filters in the AO camera to roughly match the wavefront sensor counts and limit differences 
in photon- noise error. 



Our reduction methods are described in ISkemer et al.l (120 10[ ). One difference is that 
in this paper, instead of having an overall "measurement error" , we often measure relative 
errors (between binary sources) separately from absolute errors (between one or both binary 
components and a spectroscopic standard). One benefit of this is that for sources with large 
(>2) mid-IR flux ratios, we can do an absolute calibration on just the brighter component 
with the best subset of our observations (i.e., close in airmass and time to bracketing spectro- 
scopic standards and in photometric conditions), and build up S/N on the fainter companion 
in other conditions, without having to frequently switch between the object and a standard. 
The method is also useful for tight binaries, where we can use a subset of the data to do the 
absolute calibration on the combined binary, and then use the full dataset to split the binary. 
We use this method for all but GG Tau, and XZ Tau (2009 dataset), which were observed 
only briefly, so that all of the data were used for both relative and absolute calibrations. 
The exposure times for our absolute calibrations (a subset of our full data) and our relative 
calibrations (our full data) are listed in Table 1. 
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We used t he mid-IR standard e Tau for all of our binaries except GG Tau, for which we 



used HR 1684 f ICohen et al.lll999l ). In all cases except GG Tau and XZ Tau (2009 dataset) 



we were able to measure telluric calibration errors (as in lSkemer et al.ll2010l ) with bracketing 
standards. In the cases of GG Tau and XZ Tau (2009 dataset), which were both taken in 
photometric conditions with well-matched single spectoscopic standards in both airmass and 
time, we conservatively assume 3% local errors and 5% global flux errors outside of telluric 
ozone. 

Most of the binaries in our survey have a large enough separation that splitting them is 
trivial with the 6.5 meter MMT and its adaptive optics system. However, a few of our binaries 
(GG Tau, GN Tau, XZ Tau, which are separated by 0.25", 0.41" and 0.29" respectively) 
are near or below the diffraction limit (A/-D =0.32" at lO/xm) of the 6.5 meter MMT. 
Because of MMTAO's high-Strehl and stability in the mid-i nfrared, it has resolv ed binaries 



as tight as 0.12" at 10/im using superresolution techniques (jSkemer et al.ll2008l ). By using 
our spectroscopic standard as a 1-D PSF, we fit a double (binary) PSF, for each wavelength 
of our spectrum image, fixing the separation of the binary based on imaging data. This gives 
us a measurement of the flux of each component relative to our spectroscopic standard. 

For our tightest object, GG Tau, we also flxed the curvature of the grism trace to that 
of the PSF star's grism trace, which is legitimate as long as the object and PSF star are 



obser ved at a similar airmass to avoid the effects of atmospheric dispersion (jSkemer et al. 



20091 ) . Both GG Tau and its PSF star were observed at 1.04 airmasses. 



For GN Tau, the A and B components are of similar brightness in the optical, and 
because of their small separation, the AO system's wavefront sensor was somewhat less 
effective in the direction of the binary. As a result, the binaries' PSF was slightly wider in 
the direction of the binary than the spectroscopic standard's PSF. We found that convolving 
our spectroscopic standard's PSF with a 0.05"-o" Gaussian greatly improved the residuals of 
the binary fit, although this smoothing did not significantly affect the measurement values 
of the fit. 

Spectrally dispersed images for each of our seven binaries are shown in Figure [TJ Each 
frame shows both components of a binary, with the primary on top and the secondary on 
the bottom (primary being defined as the brighter component in the optical). A vertical cut 
through each image is shown to the right. Our tightest binaries, GG Tau Aa-Ab and XZ Tau 
A-B are separated by less than the MMT's ~0.32" diffraction limit. However, the vertical 
cut through the spectrally dispersed images shows a wider spectrum, and with knowledge of 
the binary separation (obtained from acquisition images) and our high-Strehl (~98%) PSF, 
we are able to spatially resolve the spectra. 
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Spectroscopy for each of our seven binaries (plus a second epoch on one of the binaries, 
XZ Tau), is presented in Figure [2l In all cases, the A component is to the left and the 
B component is to the right. The error bars that coincide with the Earth's 9.7 fim ozone 
feature reflect the larger calibration uncertainties of our data at these wavelengths. 

A summary of all (our survey plus literature) spatially and N-band spectrally resolved 
binaries/multiples in the Taurus- Auriga star forming region is presented in Tabl e 2. We rebin 
publi shed spectra from other instruments (DP T au A-B and FV Tau A-B from lHonda et al. 
20061 and GI-GK Tau from IWatson et al.l |2009| ) to the wavelength sampling used in our 
reductions in orde r to keep our results consistent. The uncertainties for DD Tau A-B and FV 
Tau A-B are from lHonda et al.l (I20061). The uncertainties for GI-GK Tau are dominated by 
our binning method (see lSkemer et al.ll2010l ) because of Spitzer^s high sensitivity. Flux loss 
due to Spitzer slit mispointings become irrelevant in our later analysis of silicate equivalent 
widths. 



3. Analysis 



3.1. Correction for Extinction 



Our mid-IR spectra must be correcte d for extinction to avoi d contamination by the 



silicate feature of the ISM extinction cu rve fJRieke fc Lebofsky 



Class II objects are extincted by Ay <3 ( JKenyon fc HartmannI 



■19851) . In Taurus, where most 
19951 ) , the effect is not as great 



as in other star forming regions. However, som e of our objects have larger extinctions (the 
la rgest of which is FV Tau AB, with Ay = 5.3; IWhite fc Ghezll200ll ). As shown in Figure 1 
of iFurlan et al.l ( l2009l ). these extinctions will have a strong influence on the resulting silicate 
emission feature if left uncorrected. 



We correct for extinction in our measured mid-IR spectra using Ay values from lFurlan et al 



( 120091 ) and references within (listed in Table |3]), which will allow us to ma ke a consistent 



comparison between our binaries and t he single stars of iFurlan et al. (120091) in Section 4.1. 
We adopt the same extinctio n laws as iFurlan et al.l (|2009[ ): iMathij ( Il990l ) with Ry = 5 
for objects with Ay < 3 and iMcClurd ( 120091 ) for Ay > 3. Since our objects are binaries. 



we correct extinction in both components equall y, as interstellar ext i nction is thought to 
be th e same between tight binary components (IWhite fc GhezI l2001t Hartigan fc Kenvon 



(I2OO9I ) 



2OO3I ). The one exception is for our widest binary, GI-GK Tau, where iFurlan et al. 
cite separate extinction values for each component. 

Figure |3] shows our extinction corrected spectra. There are several examples of binaries 
that have strikingly similar silicate features, both in terms of their strengths and their shapes. 
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Notable example include FX Tau, GN Tau, IT Tau, DD Tau, FV Tau and GI-GK Tau. We 
quantitatively compare the strength of silicate features in binary pairs, using equivalent width 
measurements, in Section H?T1 More detailed comparisons of their detailed silicate properties 
would be useful for determining whether individual dust components (for example, certain 
crystalline species) are correlated between binary pairs, but that test is not addressed further 
in this paper. 

Our sample shows a diversity of silicate properties. Strong features (implying small 
dust grains) are present in GI Tau, GK Tau, RW Aur B and GG Tau Ab, among others. 
Weak features (no small grains) appear in IT Tau A and B and most notably in XZ Tau 
B, which shows negligible emission in 2009. XZ Tau, which was observed twice, shows 



strong variability, as was seen in multi-epoch Spitzer spectroscopy (iBary et al.ll2009l ). Our 
data indicate that silicate variability is present in both sources, and because of the strong 
variation we observe in the flux ratio of the binary, both silicate features likely affect the 
unresolved Spitzer variability. 



3.2. Equivalent Width Measurements: A Signature of Grain Growth 

Numerous authors have noted the correlation between the shape and the strength of 
the 10 micron silicate feature, where a strong feature generally implies the presence of sma ll 



grains. Since the ISM is composed of small, amorphous silicate grains ( iKemper et al.ll2004l ). 
it is assumed that circumstellar disks where these grains are observed have not undergone 
dust processing/grain growth, while disks with larger and/or more crystalline grains have 
begun grain growth. The 10 /xm silicate feature is only sensitive to a subsample of dust 
in the circumstellar disk: silicate dust grains between ~0.1/xm and ~5/im, that are at the 
optic ally thin surface of the disk and at the proper temperature to be in the silicate emission 



zone (iKessler-Silacci et al.l 120071 ). While the lO/im silicate feature is being used to probe 
grain properties in the upper-layers of the circumstellar disk, millimeter wavelengths can 
probe the cooler and deeper layers of the disk, which contain larger mm-cm sized grains. 
The shape of the millimeter continuum slope can be used to probe a characteristic mid- 
plane grain size , whic h has also been used to infer grain-growth. In a highly important work. 



Lommen et al.l (120101 ) demonstrated a correlation between the strength of the silicate feature 
(using equivalent width) and the slope of the millimeter continuum. This means that despite 
the limited emission zone of the lO/xm silicate feature, it does appear to be a marker for dust 
grain-growth throughout young circumstellar disks. This result has practical importance, 
because there are many more measurements of lO/xm silicate features than millimeter slopes, 
due to the decreasing brightness of these objects with wavelength. Binarity, in particular, is 



best studied at shorter wavelengths where the system can be spatially resolved. 



Lommen et al.l ( 120101 ) 's result correlates the equivalent width of the silicate feature with 
dust grain sizes iii the disk midplane. Similarly, silicate de-composition tec hniques (see 
Sargent et al.ll2009l : IJuhasz et al.ll2009l : lOlofsson et al.ll2010l : lOliveira et al.ll201ll and a review 



Watson 



200 



9|) establish that the silicate feature equivalent width is connected to the size 



by 

of dust grains in the upper layers of the disk. Based on these two correlations, we use 
the equivalent width of the silicate feature as a proxy for grain-growth for the rest of this 
paper. While silicate feature de-composition techniques still contain more information than 
this simple approach, equivalent width is a robust, single- number statistic, that is easily 
measured with ground-based systems that are limited in sensitivity, spectral range, and 
spectral resolution. 

We calculate equivalent width by integrating (F - Fcontinuum) / Fcontinuum, where the 
continuum is a linear fit of our flux measurements between 8.07-8. 27Atm and 12.56-12. 95/im. 
Errors are calculated using the Monte Carlo approach described in lSkemer et al.l (J2010l ). The 
dominant error source for most of our MMTA0/BLINC-MIRAC4 data is the spectroscopy at 
the edges of the silicate band, which define the (linear-fit) silicate continuum. These regions 
are particularly faint from the ground because of the low atmospheric transmission < 8/xm 
and > 12.5/im combined with decreasing flux at the longer wavelengths. The other difficult 
region in the N-band is the 9.7/im telluric ozone feature, where our absolute calibrations 
tend to be much worse than elsewhere in the spectrum (see Figure [2]). For the purposes of 
calculating equivalent width, we interpolate over the ozone region. 

The results of our equivalent width measurements of binaries, as well as 8.1 and 12.7 /xm 
continuum photometry and colors are presented in Ta ble [3l We include our MM TAO /BLINC 



MIRA C4 measurements, as well as the binaries from iHonda et al.l (120061 ) and IWatson et al. 



(I20091 ). which we interpolate in the same fashion as the MMTA0/BLINC-MIRAC4 for con- 
sistency. We do not include T Tau, UY Aur or XZ Tau, which are not used in Section 14.11 
because they are self-extincted and variable. All of the quantities listed are calculated with 
the extinction corrected spectra. 

Our ground-based equivalent width measurements use a smaller than desirable wave- 
length range (~8-12.5/im) due to truncation by atmospheric transmission. Space-based 
equivalent width measurements, moti vated by the true wi dth of the silicate feature, use 
larger wavelength ranges (for example. iLommen et al.l (120101 ) use 7.5-13. 0/im). Additionally, 
space-based equivalent width measurements often use polynomial continuum fits instead of 
linear fits and do not interpolate over telluric ozone, as is done for our ground-based measure- 
ments. For these reasons, it is worth checking that the equivalent widths calculated using 
our ground-based prescription are well-correlated with the equivalent widths calculated us- 
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ing typical space-based prescriptions. We use a set of single stars in Taurus (described in 
Sectio n 14.11) . whose extinction corrected equivalent widths are published by iFurlan et al. 
( 120091 ) based on Spitzer /lYi'^ spectroscopy. We obtained the Spitzer data, and applied our 
MIRAC4 binning procedure and equivalent width measuremeri t s pres cription. A comparison 
of our measurements with the measurements of iFurlan et al.l (120091 ) is shown in Figure HI 
The space-based prescription (true) equivalent widths are well-correlated with the ground- 
based prescription (truncated) equivalent widt hs, with the excepti on of one outlier (DM Tau, 
which has a large quadratic continuum term (jWatson et al.ll2009l ) that is not accounted for 
in the ground-based prescription). This implies that ground-based (truncated) equivalent 
widths are reliable proxies for the space-based (true) equivalent widths. 



4. Statistical Tests 



4.1. Is Grain Growth in Binaries Correlated? 



Properties shared between binary stars, such as age, composition and formation history, 
might be linked to the growth of dust grains in young stars. We test for this by determining 
if the silicate features of binary stars are more similar than the silicate features of randomly 
selected pairs of single stars. We use the ratio of extinction corrected equivalent width 
measurements (Table |3]) to quantify the similarity of two silicate features, using the reciprocal 
of the ratio when the ratio is < 1. 

The binary stars are drawn from our sample of spatially and spectrally resolved N- 
band spectra in Taurus- Auriga (Table [2]), where we exclude systems that contain infrared 
companions, that are thought to be self-extincted and variable (T Tau, U Y Aur, XZ Tau, 
as shown by Ivan Boekel et al.l l2010l : ISkemer et al.l l2010l : iBary et al.l |2009|) . This leaves us 
with 9 binaries for the test. For the pairs of single stars, we use the e x tincti on corrected 
equivalent widths of 26 Taurus-Aurigae sin gle stars from iFurlan et al.l ( 120091 ). where the 
single stars are listed in iFurlan et al.l ( 120061 ). We exclude sources with unknown and high 
[Ay > 8) extinctions, and we avoid DG Ta u, which is known to have a variable silicate 



feature (Wooden et al 



2OOOI : iBary et al.ll2009l ) and HK Tau, which has a self-extincting disk 



( IStapelfeldt et al.lll998l ). We recalculate the extinction-corrected equivalent widths of the 
single stars using our ground-based prescription (see Section 13. 2p and list these values in 
Table HI The single stars and binaries are similarly distributed on an HR diagram (Figure 
|5]), allowing us to make a reasonably unbiased comparison. 

We do an N=10'' Monte Carlo simulation to randomly pair single stars and to randomly 
sample the error bars for our binary equivalent width ratios (a histogram of the simulation 
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output is shown in Figure [6l In 90% of our trials, the average (mean) of the equivalent width 
ratios of the 9 binaries is closer to 1 than the average equivalent width ratios of 9 randomly 
selected pairs of single stars. Thus, with 90% confidence, we find that silicate features in 
binary pairs are more similar than silicate features in randomly paired single stars, implying 
that grain growth is correlated in binaries. 

There is one systematic effect that is likely to have increased the observed differences 
between the silicate features of the single star pairs. The extinction me asurements for eac h 
single star correspond to an approximately 15% error in equivalent width (JFurlan et al.ll2009l ). 
This will cause a residual differential extinction between the single star pairs, that should on 
average make the single star silicate features look more different than they actually are. We 
add 15% more uncertainty to the single star equivalent widths in our Monte Carlo simulation 
and find that our confidence is only changed by 2%. Thus, the effect of residual extinction in 
the single star pairs is small. Similiarly the extinction corrections for the binaries, described 
in Section 3.1, assumes that the binaries share the same interstellar extinction. If this 
assumption is incorrect there would be a residual differential extinction between the binary 
pairs, that should on average make the binary silicate features look more different than they 
actually are. This bias is similarly small to the single star bias. 

Our single star sample includes 3 transitions disks (GM Aur, DM Tau a nd LkCa 15), 



which have larger than normal silicate equivalent widths (see Figure H] and iFurlan et al. 



20091 ). We include them in our single star sample because some of our binaries have similarly 
large silicate equivalent widths (for example GG Tau Ab, which has a "truncated" EqW = 
3.33 ± 0.40). Excluding the transition disks (while including GG Tau) would change our 
confidence that grain growth is correlated in binaries to 83%. 

Of our 9 binaries, GG Tau Aa-Ab has the largest disparity between its silicate fea- 
ture equivalent widths. GG Tau A is a 0.25" binary with circumstellar disks around each 
compone nt and a spatially r esolved circumbinary disk (JDutrey et al.lll994l ) that is unusually 
massive (jjensen et al.l Il996l). Some authors have reported tentative detections of a dust 
streamer onto GG Tau Ab fJRoddier et al.lll996l : iPietu et al.ll201ll). which might be r eplen- 



ishing the small dust grains detected in our silicate spectroscopy (jSkemer et al.ll2010l ). The 
prese nce of streamers ha s also been inferred in tighter binary systems via periodic accre- 
tion (JMathieu et al.lll997l : I Jensen et al.ll2007l ). If the dust steamer(s) replace the dust grains 
around each components at different rates, the silicate features would not be "coeval" as we 
have assumed throughout this analysis. Excluding GG Tau (and the transition disks) from 
our Monte Carlo simulation would change our confidence that grain growth is correlated in 
binaries to 95%. 
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4.2. Does the Coevality of Binaries Reveal a Correlation of Grain Growth and 

other Properties? 

Numerous authors have used single stars (and unresolved binaries) to determine if stellar 
and disk properties might correlate with grain-growth. The strongest correlation is betw een 
X-ray luminosity and dust grain crystallinity and size ( iGlauser et al.ll2009l : iRiazI 120091 ) . A 
weak er and possibly dependent correlation might exis t between stellar mass and dust-grain 
size (JKessler-Silacci et al.l l2006t iPascucci et al.l 120091 ) . It is possible that age and/or for- 
mation history effects are obscuring the correlations between grain growth and stellar/disk 
properties. Here, we use the coevality of binaries to effectively remove the age/formation 
history affects and test correlations between grain-growth and various other parameters. 



As in Section 14. H we use the 9 binaries from our sample that do not have edge-on 
disks/self-extinction. For each of these binaries, we have compiled a table of spatially resolved 
ancillary stellar properties (Table E]) and a table of infrared flux ratios for determining disk 
colors (Table E]). We plot 6 different stellar/disk properties versus binary equivalent width 
ratios in Figure [71 With our small sample size of binaries (9), doing too many correlation 
tests would result in spurious correlation detections. Thus, the properties we have chosen 
are intended to have the most plausible links to grain-growth. Spectral type is an obvious 
catch-all, that for a fixed age, is related to the temperature, luminosity and mass of the star. 
Bolometric magnitude relat es to the total flux hitting the dust-grains, and should define 
the "silicate emission zone" ( iKessler-Silacci et al.l 120071 ). Stellar mass relates to the gravity- 
induced sinking of large dust grains to the disk midplane and disk dynamics in general. 
Ha relates to accretion, which will produce UV photons that can anneal the dust. K-L 
colors relate to the geometry of the disk's puffed-up inner-rim. K-N colors relate to the disk 
geometry at the location of the silicate emission zone. 

For each frame in Figure [TJ the x-axis shows a difference/ratio measurement for a 
star/disk property. For example, frame (a) shows the difference in spectral type between the 
A component and the B component of the binary. In all frames, the y-axis shows the ratio 
of the equivalent width ratios. Vertical dotted lines show equal spectral types, luminosities, 
masses. Ha equivalent widths, and colors. Horizontal dotted lines show equal silicate equiv- 
alent widths. This divides each frame into quadrants. We do a coarse correlation analysis 
by determining how many binaries are in the top-right or bottom-left quadrants versus the 
top-left or bottom-right quadrants. This allows us to easily answer questions of the form "Do 
the stars with earlier spectral types have larger or smaller silicate feature equivalent width?" 
Since we are dealing with coeval binary pairs, this directly relates to whether spectral type 
is correlated with grain-growth. 



In 7 out of 9 binaries, the star with the earlier spectral type has a smaller silicate equiv- 
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alent width than the star with the later spectral type. Using the binomial distribution with 
p = 50%, the probability that 7 or more out of 9 binaries would have larger equivalent widths 
in either the earlier-type components or the later-type components is 18%. This probability 
is not low enough to rule out the null hypothesis that spectral type is not correlated with 
silicate feature equivalent widths (especially given that it is 1 of 6 tests shown in Figure 
|7]). No other parameter in Figure [7] has a stronger correlation signal than spectral type. 
Additionally, none of the parameters in Figure [7] have an obvious visual trend. Thus, we do 
not find evidence for any correlations with equivalent width, other than the marginal one 
observed for spectral type. 



5. Discussion and Conclusions 

Our mid-infrared adaptive optics survey has produced 8 new spatially resolved sili- 
cate spectra of you ng Taurus-Aurigae binaries (including UY Aur, which was published in 



Skemer et al.l 120101 ). This effectively triples the number of spatially and spectrally resolved 
binaries in the Taurus- Auriga cluster, and provides the first opportunity to draw statistical 
conclusions about the spatially resolved silicate features of binaries in any young cluster. 
Our fundamental questions are (1) Is grain growth correlated between binary pairs? and (2) 
Does removing this effect reveal correlations between grain growth and other properties? 

We find with 90% confidence (when including GG Tau) or 95% confidence (when exclud- 
ing GG Tau) that the silicate features of binaries are more similar than the silicate features 
of randomly paired single stars (our ambivalence in including GG Tau stems from its large 
circumbinary disk, which might be replenishing the dust in the inner disks via streamers). 
The similarity of the silicate features in our binary pairs implies that one or more shared 
binary properties (such as age, composition or formation history) plays an important role in 
dust grain evolution. 

One possible explanation for our findings is that the environment in which a star forms 
affects its potential to grow dust grains, and eventually planets. It is a fundamental predic- 
tion of core-accretion that giant planets form more efficiently in metal -rich systems, ba. sed 



on the increased solid material available for building planetary cores ( llda fc Linll2004al Jbl). 
However, the core-accretion scenario is predicated by small dust grains coagulating quickly 
enough that aerodynamic drag does not remove too many meter-sized obje cts from the sys 



tem (a problem, which is commonly known as the "meter-sized barrier"; IWeidenschilling 



19771 ). iBrauer et al.l (|2008[ ) have shown that the removal of meter-sized objects due to radial 
drift can be overcome by increasing the disk's initial dust-to-gas ratio from 1% to 2%. The 
result is a natural consequence of increasing the rate of grain growth by decreasing the time 
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between collisions. Since binary stars form from the same fragmented core material, it is 
possible that grain growth in binary stars will proceed at similar rates in each component 
due to their shared initial gas-to-dust ratio. 

Binarity itself might also affect the silicate features of young stars. Circumstellar disks 
can be dynamically perturbed by a cor npanion star, which introdu ces spiral structure and 
truncates the radial extent of the disk (jArtymowicz fc Lubowlll994i ). In this configuration, 
l arge particles collide and fragment more frequently, which speeds up dust gra in evolution 



(IDubruUe et al.l 119951 ) while decreasing the maximum particle size in the disk (jZsom et al. 



2OIII ). Dust streamers from a wider cir cumbinary disk can also repl enish solid material in 
one or both of the circumstellar disks (jArtymowicz &: Lubowlll996l ). The dust streamers 
are likely to con tain un-evolved, ISM -like dust grains that dominate the appearance of the 
silicate feature (jSkemer et al.l I2OIOI ) . There are two young Ta urus- Aurigae binaries with 



direc tly imaged circumbinary disks: GG Tau ( jDutrey et al.lll994J ) and UY Aur (JDutrey et al. 



19961 ). In each case, we find that one of the compo nents has an un u sually prominent silicate 
feature, indicating the presence of small grains. iPascucci et al.l (120081 ) were not able to 
detect a difference between the sihcate features of (unresolved) medium-separation binaries 
and single stars, which suggests that any effect must be subtle or uncommon. However, 
spatially resolved studies of the silicate features in young binaries might reveal different 
dust-grain growth mechanisms. 

In Section 14. 2^ we found that seven of our nine binaries (i.e. 82% confidence of a 
trend using the binomial distribution) have larger silicate equivalent widths (smaller grains) 
in the secondary than the primary. This pattern was also observed in two binaries by 
Przygodda et al.l (120031 ). Assuming this trend is confirmed by larger samples of spatially 
resolved silicate spectroscopy, the cause could be based on the dynamics of binaries or simply 
a spectral-type/mass/luminosity correlation with the speed of dust grain growth. The former 
explanation could be tested by seeing if the correlation is connected with binary separation. 

Our study has been limited by its small sample size, due to the complexity and limi- 
tations of ground-based 10 fim. spectroscopy and adaptive optics, which were necessary to 
build up the current sample. Larger ground-based telescopes, such as the LBT and ELTs 
along with large space-based telescopes, such as JWST will be able to dramatically increase 
the sample of spatially resolved silicate spectra of young binaries. 
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Table 1. MMTA0/BLINC-MIRAC4 Observations of T Tauri Binaries 



Binary 


Dates (UT) 


Relative Calibrat 


ion 


Absolute Calibration 






On- 


-source Time 


{sr 


On- 


-source Time (s)'^ 


DKTau 


2010 Jan 3 




2190 






1260 


FX Tau 


2010 Jan 2 and 6 




1200 






570 


GG Tau A 


2009 Oct 2 




1200 






1200 


GN Tau 


2010 Jan 2 




2550 






750 


IT Tau 


2010 Jan 3 and 4 




2385 






2385 


RW Aur 


2010 Jan 5 and 6 




765 






630 


XZ Tau 


2009 Jan 14 




160 






160 


XZ Tau 


2010 Jan 5 




1230 






270 



''Relative fluxes between binary pairs were measured using all available data. Abso- 
lute calibrations (using the brighter star or the "un-resolved" binary pair) were mea- 
sured using the best subset of data to minimize errors caused by changing atmospheric 
conditions. 
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Table 2. N-band Spatially and Spectrally Resolved Taurus Binaries 



Binary 


Separation 


(") 


Separation Ref. 


N-band spectra Ref. 


DD Tau A-B 


0.56 




1 


2 


DK Tau A-B 


2.32 




1 


this work 


FV Tau A-B'' 


0.71 




1 


2 


FX Tau A-B 


0.89 




1 


this work 


GG Tau Aa-Ab'' 


0.25 




this work 


this w^ork 


GI-GK Tau<^ 


12.9 




3 


4 


GN Tau A-B 


0.41 




this work 


this work 


IT Tau A-B 


2.4 




5 


this work 


RW Aur A-Bd 


1.40 




1 


this work 


T Tau N-Sa-Sb 


0.69/0.13 


6 


7 


UY Aur A-B 


0.88 




1 


8 


XZ Tau A-B° 


0.29 




this work 


this work 



Note. — A hst of T Tauri binaries in the Taurus- Auriga star-forming regions with 
spatially resolved N-band spectroscopy. We onl y include binaries separated by < 30", 
which are then likely to be bound and coeval JKraus fc Hillenb rand 20091). Due to 
non-negligible orbital motion in the cases of GG Tau, GN Tau and XZ Tau, astrom- 
etry is measured from 8.7/im acquisition images preceding our N-band spectroscopy. 
We use a Gemini as our astrometric standard. T Tau Sa-Sb also h as fast orbital 
motion, which is parameterized, most recently, bv lKohler et al.l 1120081 ) 

"■does not include FV Tau /c, which is a 0.7" binary located 12.3" from FV Tau 
jWhite fc GheJboOll) 

^does not include GG Tau B, which is a 0.25" binary located 10.3" from GG Tau 
A l lWhite fc GhejboOlh 

'^GK Tau is itself a 2.4" binary JHartigan et ahlllQM) . However, GK Tau B was 
not detected in our 8.7 micron imaging. 

dRW Aur was found to be a triple by JGhez et al.ll993h . but llWhite fc GheJbOOlh 
claim the 3rd star was probably a false detection. 

"^X Z Tau A is itself a 0.09" binary at mm wavelengths lICarrasco-Gonzalez et al.l 
l2009h . 



References. — (1) IWhite fc Ghe3 ll200l[): (2) iHonda et al.l ll2006h : (3) 
iHar tigan ct~all lll994l); (4 ) iFurlan et al.l ll2q06h. l\yatson et al.l ||2009|) and o thers; 
(5) iDuchene et al.l lll999h : (6) iKohler et al.l ll2008h ; (7) iRatzka et al.l ll2009h ; (8) 
ISkemer et al.l ||2010|') 
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Table 3. Extinction Corrected Silicate/Continuum Measurements of T Tauri Binaries 





Avimag)"^ 


Silicate EqW (Aim) 


8.1 firn^ flux (Jy) 


12.7 fim'= flux (Jy) 


DD Tau A 


1.0 


0.95±0.06 


0.45±0.01 


0.73±0.00 


DD Tau B 


1.0 


1.53±0.11 


0.34±0.01 


0.43±0.01 


DD Tau A/B 




0.62±0.06 


1.33±0.06 


1.72±0.03 


DK Tau A 


1.3 


2.80±0.18 


0.66±0.03 


0.78±0.02 


DK Tau B 


1.3 


1.76±0.25 


0.15±0.01 


0.14±0.01 


DK Tau A/B 




1.60±0.18 


4.51±0.14 


5.52±0.30 


FV Tau A 


5.3 


1.03±0.10 


O.SOitO.OO 


0.32±0.01 


FV Tau B 


5.3 


1.44±0.11 


0.43±0.01 


0.54±0.01 


FV Tau A/B 




0.72±0.09 


0.70±0.03 


0.59±0.02 


FX Tau A 


2.0 


2.61±0.39 


0.14±0.01 


0.17±0.02 


FX Tau B 


2.0 


2.34±0.43 


0.05±0.00 


0.07±0.01 


FX Tau A/B 




1.12±0.26 


2.91±0.22 


2.45±0.35 


GG Tau Aa 


1.0 


0.91±0.25 


0.36±0.02 


0.39±0.04 


GG Tau Ab 


1.0 


3.33±0.40 


0.24±0.02 


0.28±0.02 


GG Tau Aa/Ab 




0.27±0.08 


1.47±0.10 


1.38±0.14 


GITau 


2.3 


1.93±0.05 


0.68±0.01 


0.79±0.00 


GKTau 


1.1 


2.76±0.09 


0.69±0.02 


0.86±0.00 


GI Tau/GK Tau 




0.70±0.03 


0.99±0.03 


0.92±0.00 


GN Tau A 


3.5 


1.04±0.21 


0.14±0.01 


0.15±0.01 


GN Tau B 


3.5 


1.66±0.20 


0.23±0.01 


0.21±0.01 


GN Tau A/B 




0.62±0.12 


0.59±0.03 


0.75±0.04 


IT Tau A 


3.8 


0.86±0.18 


0.20±0.01 


0.17±0.01 


IT Tau B 


3.8 


1.06±0.37 


0.08±0.00 


0.07±0.01 


IT Tau A/B 




0.81±0.26 


2.51±0.07 


2.29±0.27 


RW Aur A 


0.5 


1.10±0.07 


1.05±0.05 


1.26±0.05 


RW Aur B 


0.5 


2.45±0.48 


0.12±0.01 


O.lliO.Ol 


RW Aur A/B 




0.45±0.10 


8.43±0.62 


11.81±1.59 



Note. — Afl error bars are calculated by taking the interquartile range of our Monte Carlo simu- 
lation, and dividing by 1.34. This is a robust estimate of a Gaussian sigma. 

'^Av values from lFurlan et al.l ||2009| ') 

''Flux measured between 8. 07^101 and 8.27Aim 

'^Flux measured between 12.56/im and 12.95/xm 
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Table 4. Single Stars Used in Comparison with Binary Stars 



Av('mag)^ Spectral Type'' Silicate EqW (Mm)'' 

MO 0.21 

K7 1.18 

K7 1.89 

K7 1.61 

K3 0.69 

MO 1.05 

K7 0.42 

MO 1.21 

K7 0.40 

Ml 1.43 

MO 0.51 

MO 0.65 

M0.5 1.06 
0.80 

K5 1.66 

M4 1.75 

M5 1.25 

K3 3.50 

MO 1.49 

K3 1.69 

MO 2.63 

M0.5 1.26 

K5 3.98 

Gl 4.02 

M6 0.99 

K7 2.02 

^Av values and spectral types from lFurlan et al.l 1120091 ) 

''Eq uivalent Widths are based on Spitzer spectra from iFurlan et al.l 
1 120091) and are calculated using our ground-based prescription as de- 
scribed in Section [XT] 



04108-1-2910 


1.4 


AA Tau 


1.8 


BP Tau 


1.0 


CI Tau 


2.0 


CWTau 


2.8 


CXTau 


1.3 


CYTau 


1.7 


DETau 


1.2 


DLTau 


1.5 


DM Tau 


0.7 


DNTau 


0.6 


DO Tau 


2.0 


DP Tau 


0.6 


DRTau 


1.2 


DS Tau 


1.1 


F04147-t-2822 


2.5 


FN Tau 


1.4 


GM Aur 


1.2 


GO Tau 


2.0 


HP Tau 


2.8 


IP Tau 


0.5 


IQ Tau 


1.4 


LkCa 15 


1.2 


RYTau 


2.2 


V410 Anon 13 


5.8 


V836 Tau 


1.1 
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Table 5. Ancillary Stellar Properties of Taurus Binaries 



Binary 


Spectral Type 


Ha, EqW (A) 


Ref. 


Mtoi (mag)'' 


Mstar/MQ'^ 


DD Tau A 
DD Tau B 




M3.5 
M3.5 


206 
635 


1 


5.08 
5.40 


oAstn 


DK Tau A 
DK Tau B 




K9 

Ml 


31 
118 


2 


3.66 
5.35 


ny-i+O.Oy 
f, pi +0.06 


FV Tau A 
FV Tau B 




K5 
K6^ 


15 
63 


3,1 






FX Tau A 
FX Tau B 




Ml 
M4 


13 
1.0 


4 


4.91 
5.97 


69+"-0^ 
U-tl^-0.03 

n 90+0. 16 


GG Tau Aa 
GG Tau Ab 




MO 

M2 


57 
16 


1,3= 


4.50 
4.87 




GITau 
GK Tau 




K6 
K7 


20 
22 


5,6 






GN Tau A 
GN Tau B 


M2.5 


(unresolved)*^ 


59 (unresolved) 


7 






IT Tau A 
IT Tau B 




K3 

M4 


21.7 
147 


4 


3.84 
5.98 


1 0+"'^ 

J^-U-o.i 

n 90+0.15 


RW Aur A 
RW Aur B 




Kl 
K5 


76 
43 


3 


3.11 
5.01 


1.4±rr 
0.86±0;1J 



^Mfj^l and Mstar are calculated bv lKraus fc Hillenbrandl 1120091 ) 

''Note that iHartiean fc KenvonI 1120031 ) found no pliotospheric absorption lines in their STIS 
spectrum of FV Tau B. 

HHartigan fc KenvonI 1120031) 's GG Tau spectr um is saturated at Ha , so we use Ha equivalent 



width measurments from lWhite fc Ghea 11200 ll 

"* Although GN Tau does not h ave a published resolved spectral typ e, it has been resolved 
in the optical llSimon et al.lll996h . near-infrared l lWhite fc Gheal200ll) and mid-infrared (this 
work). Throughout this range of wavelengths, it has a near equal flux ratio, so it is likely that 
the binary consists of two nearly equal mass stars. 



Ref erenc es. — (1 
(I2OOII'): (4' )lDuchene et al 
(Tl lWhite fc BasrU.200: 



iHartigan fc KenvonI IJOOJ): (2) iMonin et al.l lll998h: (3) IWhite fc Ghe3 



lean & KenvonI i:^00a): VI) IMonm et al.l 1199!^): l3)IWhite&: Ghe2l 
il999h ; (5) JKenvon fc HartmannI lll995h - (6) JHartigan et al.l 1I199I) : 
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Table 6. Flux Ratios of Taurus Binaries 



Binary 


K Flux Ratio 


L Flux Ratio 


Near-IR Ref. 


12.7 fim Flux Ratio 


DD Tau A-B 


1.17±0.06 


1.72±0.13 


1 


1.72±0.03 


DK Tau A-B 


4.21±0.30 


4.02±0.13 


1 


5.52±0.30 


FV Tau A-B 


1.49±0.01 


0.69±0.02 


1 


0.59±0.02 


FX Tau A-B 


2.20±0.02 


2.97±0.03 


1 


2.45±0.35 


GG Tau Aa-Ab 


1.95±0.16 


2.23±0.01 


1 


1.38±0.14 


GI-GK Tau'' 


0.77±0.03 


0.77±0.07 


2,3 


0.92±0.00 


GN Tau A-B 


1.17±0.03 


1.25±0.03 


1 


0.75±0.04 


IT Tau A-B 


6.37±0.76 


16.0±4.4 


1 


2.29±0.27 


RW Aur A-B 


4.23±0.16 


6.16±0.24 


1 


11.81±1.59 



''Photometry for GI Tau and GK Tau have been extinction corrected individually using the 
Ay values listed in Table [3] 



Ref erences. — (1) IWhite fc Ghe3 ll200ll '): (2) ISkrutskie et alj 1I2OO6I '): (3) iLuhman et"al] 

Hooi) 
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Sep = 0.25' 




Fig. 1. — Spectrally dispersed images of 7 T Tauri binaries that we have spatially and 
spectrally resolved for the first time at N-band. The wavelength range (shown at the top) 
is 7.9/im to 12.5/xm and has a suppressed flux at 9.7/im due to telluric ozone. The images 
are each 4.4" in the spatial direction and are displayed with (different) linear stretches. At 
the right of every spectrally dispersed image is a vertical (spatial direction) profile cut. GG 
Tau and XZ Tau are separated by less than the diffraction-limit of the 6.5 meter MMT but 
are superresolved due to the excellent stability of the MMT's adaptive optics PSF in the 
mid-infrared. This is evident in the vertical profiles, which show that that GG Tau and XZ 
Tau are wider than the PSF (seen in the well-separated binaries below). 
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Fig. 2. — Spectra of 7 T Tauri binaries that we have spatially and spectrally resolved for 
the first time at N-band. XZ Tau is shown twice because we observed it once in 2009 and 
once in 2010. In all cases, the primary is shown on the left and the secondary is shown on 
the right. The grey region on each plot shows the la range of our continuum fit, which we 
use when calculating equivalent widths. Continued 
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Fig. 2. — Continued 



-24- 



1.5 




' 


o >. 


.^#^ 


^^ 


Cx 1.0 






■>< = 




"^^ 


DK Tou 
correc 

o 






0.0 




_ 



9 10 11 12 
Wovelength (/tm) 



13 




9 10 11 12 
Wovelength (jim) 




9 10 11 12 
Wavelength (/im) 



0.20 [ 



o >v 

t^ 0.15- 



Itj 0.10- 

3 O 
3 V 

- t 0.05 - 
< o 
_ " 

0.00 L 



9 10 11 12 
Wavelength (/inn) 



13 



•^ X 

X 3 
0)^ 






O.b 


.,,,,,,_ 


0.5 


: ^*H 1 


0.4 


0.3 


0.2 


■; 


0.1 


r ] 


no 


■ 



9 10 11 12 
Wovelength (/im) 



13 14 




9 10 11 12 
Wovelength (/xnn) 




9 10 11 12 
Wovelength {fim) 



3 extinction 
d flux (Jy) 

o o o 


- / 


/^^ 




.^ 0.2 
oj 

zo 0.1 
0.0 






^ 1 



9 10 11 12 
Wovelength {fim) 



13 



Fig. 3. — Extinct ion c orrected spec t ra fro m Figure [2] with 3 additional binaries from 
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Fig. 4. — We take space-based equivalent width measurements from ( JFurlan et al.ll2009l ) and 
compare them with equivalent width using the same data but a different (i.e., ground-based) 
measurement prescription. The different prescriptions produce well-correlated equivalent 
widths (as shown by the best-fit, dotted line), where the ground-based (truncated) equiva- 
lent widths are systematically lower by 36.6%. This means that ground-based (truncated) 
equivalent width measurements can be used as a proxy for space-based (true) equivalent- 
width measurements. Note that our plot has one outlier, DM Tau, due to the large quadratic 
term in its continuum that is not included in our ground-based prescription. The 3 transition 
disk objects all have unusually large equivalent widths. 
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Fig. 7. — Correlations between binary star silicate feature equivalent width ratios and a set 
of ancillary stellar properties and disk colors. The dotted lines show equal properties between 
binaries. No obvious correlations are present. Frame (a) shows that in 7 out of 9 binaries, 
the secondary has a larger silicate feature equivalent width {EqWio^A/ EqWio^B < !)• The 
probability that 7 or more out of 9 binaries would have larger equivalent widths in either 
the earlier-type components or the later-type components is 18%, which is low enough to 
warrant further study, but not low enough to prove a statistical difference. 
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